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Abstract. We determine the two loop correction to Gribov's mass gap equation for quantum 
chromodynamics in the Landau gauge in the MS scheme by computing the two loop correction 
to the horizon condition derived from Zwanziger's local renormalizable Lagrangian which in- 
corporates the Gribov parameter. We verify that with the explicit result, the two loop ghost 
propagator is enchanced in the infrared. 
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In 1978 Gribov demonstrated that globally fixing a gauge in a non-abelian gauge theory could 
not be performed in a unique way, p. This is due to the fact that different gauge configurations 
could satisfy the same gauge fixing condition. Related to this was the observation that the first 
zero mode of the Faddeev-Popov operator defined the edge or horizon of a region, known as the 
Gribov volume, inside which one had to restrict the region of integration of the path integral 
defining the quantum field theory, ^ . One consequence of the existence of this restricted domain 
was that the form of the gluon and ghost propagators in the infrared regime differed significantly 
from that usually used in perturbative, and therefore ultraviolet, calculations. More concretely, 
the ghost propagator was enhanced and the gluon propagator suppressed as the momentum 
decreases. Such behaviour is believed to be an element of the confinement mechanism which 
is nonperturbative. Another feature of Gribov's original work was the definition of the volume 
of the Gribov region*. This was quantified in terms of a dimensionful parameter, known as 
the Gribov volume or mass, which was not independent but satisfied a gap equation. Gribov 
computed the explicit form of the gap equation to one loop and related it to the running coupling 
constant, pp. Subsequently, Zwanziger examined the Gribov path integral formulation of the 
problem and managed to construct a localized renormalizable Lagrangian, which included the 
Gribov mass explicitly, in addition to several new ghost fields over and above the usual Faddeev- 
Popov ghosts originating from the standard gauge fixing procedure, [SI IE]- When these 
new ghosts are eliminated by their equations of motion one recovers the non-local Lagrangian 
introduced by Gribov. The renormalizability properties of the Zwanziger Lagrangian have been 
studied in detail in [3 |Hj in the Landau gauge and there is an interesting structure. First, the 
wave function renormalization constants of the extra Zwanziger ghosts are identical to those of 
the Faddeev-Popov ghosts in the Landau gauge despite being spin-1 fields with a different colour 
structure. Second, the renormalization of the Gribov mass is not independent. Specifically the 
anomalous dimension of the Gribov mass parameter is given by a combination of the /3-function 
and gluon anomalous dimension. These have been demonstrated to all orders in perturbation 
theory using the algebraic renormalization technique, |Z1IH]- That Gribov's original formulation 
of the limitations of the quantization of a non-abelian gauge theory can be recast in a localized 
renormalizable Lagrangian means that it ought to be possible to use the Zwanziger Lagrangian, 
|S1 m El El ) to study problems in the presence of the Gribov mass. For example, it should be 
possible as an initial exercise in this direction to compute the two loop correction to Gribov's 
mass gap equation and therefore refine the relation between the Gribov mass and the coupling 
constant. This is the main aim of this article where we will concentrate on establishing the two 
loop correction to Gribov's mass gap equation in the Landau gauge in quantum chromodynamics, 
QCD. 

First, we recall the Zwanziger Lagrangian in the conventions and notation we will use. For 
QCD we have in d spacetime dimensions, jSHHElini, 
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where 7 is the Gribov mass parameter and 

= - \gI^G''^"' - —id^^A'^f - (^d'^D.c^ + iijj'^m'^ (2) 

where is the gauge field, and are the Faddeev-Popov ghosts and it is understood that 
we will only consider the Landau gauge given by a = 0. Further, the Zwanziger ghost fields, 

*For the reader interested in more detailed background to the Gribov problem, the lectures of give detailed 
working of the calculations of ^ 
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{4>1l!' , 4'1l!'} , are complex conjugate commuting fields and {uj'jf , Co'j^^} are complex conjugate anti- 
commuting fields. The remaining conventions are that = — d,yA'j^ — gf°'^'^A^^A1, g is 

the coupling constant, T'^ are the colour group generators whose structure constants are f°'^'^, 
ip^^ is the (massless) quark field, 1 < a < Na, 1 < / < Np and 1 < i < Nf with Np and Na the 
dimensions of the fundamental and adjoint representations respectively, and Nf is the number 
of quark flavours. The covariant derivatives are defined by 

{D.cp^r' = d,<pf - gr'^A^ckf . (3) 

It is worth noting the notational contrasts between this version of the Lagrangian and that of 
other articles. We have defined the Gribov mass parameter, 7, to be of mass dimension one. 
Some authors use 7^ or 7^^, in our notation, as the equivalent Gribov parameter with respective 
mass dimension two and four. Second, we do not include a coupling constant with 7 in the 
mixed mass term for the same reasons one does not include it for, say, a quark mass term. 
Third, the appearance of the factor in the mixed quadratic term plays a key role not only in 
the derivation of the propagators of but in verifying the correctness of the gap equation we 
will determine. With this version of the Gribov-Zwanziger Lagrangian, we follow |6j in noting 
that the Gribov horizon condition, which leads to the mass gap, is 

r'^{A'^^{x)<l^l^{x)) = (4) 

where all quantities are bare and @ is derived from the equation of motion of 0^'' in . 

At one loop it is elementary to see that this condition leads to Gribov's original gap equation 
of pp. For instance, from the terms quadratic in the fields of the gluon and commuting ghost 
propagators are given by 



_ siacTbd fabe fcde^i 

Xacxbd 

{i^f{p)Cjl\-p)) = -'-^v,u (5) 



in momentum space where 



p2 



PAp) = Vt^u - ^ (6) 



and f°-'^'^f^<^'^ = Ca^""^ ■ Concerning conventions here the appearance of Ca in the denominator 
factors is a direct consequence of our choice of the coefficient of the mixed quadratic term of 

In 1^ the convention was to define the mixed quadratic term to have an additional factor 
proportional to ^/Ca whence the corresponding term of © would be merely 7^. However, the 
key point in the derivation of Q is the role played by the factor of in the mixed mass term of 

Ordinarily in deriving the propagators from a Lagrangian one isolates the terms quadratic 
in the fields in momentum space and inverts the associated operator or matrix. In the case of 
this procedure is followed but not only is one dealing with a mixing between the gluon and 
commuting ghost fields but the former is real whilst the latter is complex. To deal with this 
consistently, one must write this matrix with respect to the basis ("^^^j 'i^'^)- In the absence of 
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the mixing term one would not need to do this as the matrix to be inverted is block diagonal. If, 
by contrast, the factor of was omitted from (^, then the common factor in the propagators 

of © would be [{p^Y + 2C ai'^] as was used in jS]. Since the original Gribov article has the factor 
[(p^)^ + Ca7^] we choose to include the normalization of explicitly in to have propagators 
which are consistent with those of T] . We will comment further on the significance of this factor 
later, except to note here that the appearance of in ((21) derives from the form of the basis 
chosen. Finally, concerning the mixed propagator of ©, we note that this mixing disappears in 
either the limit as 7 ^ or in the abelian limit which is formally defined as /"'"^ — > implying 
Ca ^ 0. 

Whilst this type of mixed propagator will complicate the problem of performing loop cal- 
culations with (^, at one loop it is straightforward to use the mixed A'^-(f)^J^ propagator itself 
to evaluate the vacuum expectation value of Q. Thoughout our loop calculations we will use 
dimensional regularization in d = 4 — 2e dimensions and subtract the divergences using the MS 
scheme. Using partial fractions and the standard one loop massive integral 



k [A;2 + m^] 



(m2)K 



-r (1 - \d) 
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where - j ^3^) 
equation in four dimensions of 



/ ^J^^ and m is a general mass argument, we reproduce the finite Gribov mass gap 
" as 
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a + 0{a^ 
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where a = ^^/(IGvr^) and is the MS renormalization scale introduced to retain a dimension- 
less coupling constant in d-dimensions and incorporates the usual numerical factor of 47re~^E 
with 7e the Euler-Mascheroni constant. In deriving this finite expression we have introduced 
renormalization constants for the bare objects appearing in @. Specifically these are defined 
by 
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where the subscript o denotes the bare quantity here. Since the main aim of this article is to 
extend (jHJ to the next order, we have made use of the symbolic manipulation language Form, 
[S], to first reproduce the original mass gap of P before considering the two loop calculation. 
In order to check that we have used a consistent set of Feynman rules in the Form programme, 
we have first explicitly renormalized to two loops in MS and verified that the Slavnov- Taylor 
identities derived in [71 |H] are correctly reproduced. However, such a two loop renormalization 
can be performed with the massless version of (Q). This means that the Mincer algorithm 
for massless 2-point functions, |lUj . was the tool used for this particular computation. For 
completeness we note that the explicit expressions for the renormalization constants for the 
quantities which do not ordinarily arise in the treatment of the usual QCD Lagrangian are 
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where T'^T'^ = Cpl and tr [T^T^ 
malization group functions as 



TpS""^ . The pole structure can be encoded in the renor- 



l<t>{a) = luj{a) 



hjAa + [mCATp^ - 95Cl] ^ + 0(a3) 



7^(a) = [ler^TV/ - 35C7a] ^ + f280CA7>iV/ - 449Cl + 192C7Fr^iV/l ^ + 0(a=^ 

48 L J 192 



where, in the Landau gauge, 



ainz^ 

da 



7^ (a) = /3(a) 



91nZ^ 
da 



7^ (a) 
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and in (|13j) we have used the fact that the renormahzation constants do not depend on 7. Clearly 
70(a) and 7t^(a) are equivalent to 7c(a) in agreement with the expectation of [71 |H] and 77(a) 
satisfies the Slavnov- Taylor identity derived in [3 |H] . 

Concerning the explicit expression we note that it agrees with that of ^ when evaluated 
in dimensional regularization in MS. In this respect it is important to note the role played by 
the explicit appearance of the dimension d in 0. In evaluating (jlj) we have not set d to be 4 
on the right hand side at the outset. Whilst this may appear to be a trivial point here, since 
the contribution from the 0(e) term of d = 4 — 2e is present in the constant term of (jHl), it will 
turn out that to have a finite two loop gap equation one must retain (i as 4 — 2e in Q prior to 
renormalization in MS. Having established ^ it is worth noting that for an abelian theory the 
equation is trivially satisfied since Na is formally zero in the original defining horizon condition, 
©. 

Zwanziger's elegant reformulation of the Gribov problem and the horizon definition, Q, 
in his Lagrangian immediately opens up the path to computing the two loop correction to 
(jH)) which we now detail. One simply has to evaluate the two loop corrections to the vacuum 
expectation value /"^'^(A"'^^^'^). Since this will involve two loop massive vacuum bubbles it 
remains to determine the set of Feynman diagrams and the procedure to evaluate them. One 
problem with the former is the complication of having a mixed A^-^^'^ propagator. To correctly 
establish the Feynman graphs we used the Qgraf package, using an adaptation provided 
by the author^^ of Consequently there are 17 two loop Feynman diagrams to determine. 
These were computed using a symbolic manipulation programme written in FORM which took 
the Qgraf output, substituted the Feynman rules for (0) and broke up the vacuum bubbles 
into a common form, which was 



h {jnijUi^jni^; a, b, c 



1 



(14) 



^We are very much indebted to Dr P. Nogueira for his elegant solution of the mixed propagator problem in 
Qgraf. 
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in addition to the product of two one loop integrals of the form of (0). Here {a,b,c} are 
strictly positive integers and the two loop vacuum bubbles have potentially three mass scales 
after one applies partial fractioning. In other words the set of masses {mf , ml} can take 
any combination of values in the set {0, ^7^, — i7'^}. Since the general form of H14|) has been 
evaluated to the finite part in d-dimensions for arbitrary mi, m2 and m^, it was a straightforward 
exercise to make the appropriate identifications in a Form routine using the results of |12l 1131 
E]- Though it is worth noting that whilst the combinations of {ml, 1712, m'^} may appear to 
lead to a complex value for the integrals, in the overall sum for P''%A''^'(|P^) the final result 
remained real which was a useful check. Further checks resided in the fact that the correct 
divergence structure resulting from combining the basic integrals (O and (|14j) in summing up 
all the contributions from all the Feynman diagrams was correctly cancelled by the available 
renormalization constants. In this respect we have followed the standard procedure of jl5j for 
renormalization in automatic calculations. In this the two loop calculation is performed for bare 
parameters before their renormalized versions are introduced which automatically introduce the 
appropriate counterterms. Consequently, one is left with the main result of this article which is 
the finite two loop correction to ((SJ, 
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where S2 = (2\/3/9)Cl2(27r/3), Cl2(x) is the Clausen function and Q{n) is the Riemann zeta 
function. For the interested reader S2 and ^/bC,{2) arise from the finite part of integrals such as 
l2{m'^ , m'^ , m'^ ] 1, 1, 1) and /2(m^, m^, —m?; 1, 1, 1) respectively. We also note that taking d as 4 
in Q would have resulted in an extra divergence which could not be cancelled. 

Having established the two loop gap equation for the Gribov mass, we now comment on 
its implication for the ghost propagator. In ^ Gribov showed that the one loop condition Q 
ensured that the one loop correction to the ghost propagator was enhanced in the infrared. More 
specifically if one writes the full ghost propagator as 



(16) 



then in the ultraviolet, the ghost propagator has the usual (perturbative) behaviour of 1/p^. 
However, in the infrared the propagator behaves as l/(p^)^ as ^ 0, 0. Gribov indicated 
that this was a property of a confining non-abelian gauge theory. Moreover, he computed the one 
loop correction to the ghost propagator to O(p^) and found that 1 + n(0) = provided that 7 
satisfied (jSJ. This observation of 1 + n(0) = 0, which is known as the Kugo-Ojima confinement 
condition, JH]) should also hold at two loop in if the Gribov-Zwanziger Lagrangian is a 
consistent formulation of the Gribov horizon condition. Therefore, we have computed the two 
loop correction to u{p^) by generating the 31 Feynman diagrams using Qgraf, expanding 
them to 0{p'^) and evaluating the resulting vacuum bubbles. It transpires that this 0{p'^) 
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correction to u{p'^) is exactly the same as the two loop part of (P3|) and therefore the Kugo- 
Ojima condition is satisfied at two loops precisely for all colour groups and massless quarks. 
Hence, the ghost propagator is enhanced at two loops in the Gribov-Zwanziger Lagrangian and 
has a behaviour in the infrared which is in qualitative agreement with other approaches. 

For instance, ghost enhancement has also been observed in explicit (non-perturbative) studies 
using the Schwinger-Dyson formalism and on the lattice. More specifically in the Schwinger- 
Dyson computations estimates for the exponent of the power-law behaviour of the gluon and 
ghost propagators in the infrared have been extracted and give values which are different from 
those of the ultraviolet form of the propagator, Also, in several lattice studies the Kugo- 

Ojima confinement condition itself has also been examined. See, for instance, [W \ \'2U \ nT \ \22\. 
where one recent lattice estimate for m(0) is — 0.83, |22j . 

It is worth commenting on one final aspect of the two loop calculation of (|15|) . If one uses the 
propagators of without the factor of in © then not only would the ghost propagator not 
be enhanced but at an earlier point of the study, the ghost 2-point function would not actually 
be finite either. This further justifies the derivation of © from and what would maybe 
initially appear as a peculiar convention in the definition of the mass term. More importantly, 
that (|15|) is exactly what is required for ghost enhancement provides a strong check on the result 
H15|) as well as ensuring that one has a consistent set of Feynman rules for non-zero 7. Next, 
in concentrating on the consistency aspect of the ghost with respect to enhancement, the gluon 
propagator is also expected to have infrared behaviour different from the usual ultraviolet form, 
^ . Whilst we have not examined the two loop corrections to the gluon propagator as it 
is technically more difficult to analyse than the ghost in the ^ limit, we do not believe 
its behaviour in the Gribov-Zwanziger Lagrangian, (^, to be inconsistent with the expectation 
that the propagator vanishes, pi ITS]. 

In conclusion, we have provided the O(a^) correction to the Gribov mass gap equation for 
QCD in the Landau gauge. Reassuringly the explicit form, (|15j) . guarantees the enhancement of 
the ghost propagator at two loops. Given the recent interest in both the Zwanziger approach to 
incorporating the Gribov problem in a localized renormalizable Lagrangian and in gauges such 
as linear covariant, j2,Sj . and maximal abelian gauge, |23) if would appear plausible that one 
could extend the one loop mass gap equations in those gauges to two loops as well. Moreover, 
given that |H1 I25| I26j also examined the Gribov problem using the local composite operator 
formalism to include the dimension two composite operator ^A^-^A^ it would be interesting to 
extend that one loop analysis to see whether the operator condenses and lowers the vacuum 
energy as it does in the case when the Gribov volume is regarded as infinite. Finally, we note 
that we believe that this is the first non-trivial loop computation performed with Zwanziger's 
Lagrangian. Given that the Gribov mass can be incorporated in calculations now, it would 
be interesting to examine the effects the presence of 7 has in phenomenological analyses. For 
instance, using QCD to examine deep inelastic scattering problems one will inevitably wish to 
extend such analyses towards the infrared. It is important that the Gribov limitations are taken 
into account because theoretical predictions may no longer be relevant. One case in point is the 
extraction of estimates for the dimension four condensate {C^^G"" ^'^) . Since the expansion of the 
gluon propagator in the presence of the Gribov mass naturally gives rise to a power correction of 
dimension four, it would seem important that the consequences of a non-zero 7 are understood 
in the corresponding underlying operator product expansion. 
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